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Abstract 

We consider the electrostatic inverse boundary value problem also 
known as electrical impedance tomography (EIT) for the case where the 
conductivity is a piecewise linear function on a domain fl C R" and we 
show that a Lipschitz stability estimate for the conductivity in terms of 
the local Dirichlet-to-Neumann map holds true. 


1 Introduction 

We consider the inverse boundary value problem (IBVP) associated with the 
elliptic equation for an electric potential, where the objective is to recover elec¬ 
trical resistivity, or conductivity, from partial data. We focus our attention on 
the stability of this inverse problem, in particular, when the conductivity is 
isotropic. We obtain a Lipschitz stability result if the conductivity is known to 
be piecewise linear on a given domain partition. 

We let n be a bounded domain in R", n > 2. In the absence of internal sources, 
the electric potential, u, satisfies the elliptic equation 

(1.1) div( 7 VM) = 0 in n, 

where the function 7 signifies the conductivity in fl; 7 is a bounded measurable 
function satisfying the ellipticity condition, 

(1.2) 0 < < 7 < A, almost everywhere in 17, 

for some positive A S R. The inverse conductivity problem consists of finding 7 
when the so-called Dirichlet-to-Neumann (DtoN) map 

(1.3) A.^ : u\on&Hi{dn)—>'fVu ■ v\on e H~i {dit) 
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is given for any weak solution u G to (11.11) . Here, denotes the unit 

outward normal to i9H. If measurements can be taken on a portion E of dfl 
only, then the relevant map is referred to as the local DtoN map. 

This inverse problem has different appearances, namely, as electrical impedance 
tomography (EIT) and direct current (DC) method or electrical resistivity to¬ 
mography (ERT) in geophysics (belonging to the class of potential field meth¬ 
ods). Although the mathematical framework of this paper is the one described 
by (|1.1I) - (I1.3I) . the application we have in mind is the determination of the re¬ 
sistivity p = 7 “^ in DC or ERT methods corresponding with the following type 
of experiment or “sounding”: a current is injected into the ground through a 
pair of electrodes at the boundary while the voltage is measured with another 
pair of electrodes. Thus the data, viewed as an operator, can be identified with 
the so-called Neumann-to-Dirichlet (NtoD) map. We note that in the mathe¬ 
matical literature the use of the DtoN map as the data is more common. The 
DtoN map is invertible on its range. Indeed the applied boundary current fluxes 
must have a vanishing average. We note that the solution is defined up to an 
additive (grounding potential) constant. Whereas it is well known that, at a 
theoretical level, the knowledge of either of the two maps is equivalent, matters 
may be more complicated when, in different applications, the physical settings 
provide different discrete and noisy samples of such maps. The NtoD map, upon 
applying it to a particular subset of currents, provides the so-called apparent 
resistivity. To be precise, the apparent resistivity is a geometrical (acquisition) 
factor multiplied by the ratio of voltage (potential difference) over current. 

The first mathematical formulation of this inverse problem is due to Calderon 
in the context of EIT in where he addressed the problem of whether it is 
possible to determine the (isotropic) conductivity from the DtoN map. To be 
precise, Calderon investigated the injectivity of the map 

Q ■ 7 ^ Qjj 

where Q-y{4') is the quadratic form associated to A.^, by linearizing the problem. 
As main contributions in this respect we mention the papers by Kohn and 
Vogelius |Koh-VolllKoh-Vo^ . Sylvester and Uhlmann |Sy-U| , and Nachman [N]. 
We wish to recall the uniqueness results of Druskin who, independently from 
Calderon, dealt directly with the geophysical setting of the problem in m-m- 
We also refer to [Bo| , [Che-I-N] and [U] for an overview of recent developments 
regarding the issues of uniqueness and reconstruction of the conductivity. 

It is well known that the IBVP of determining the conductivity 7 from the 
DtoN map is ill-posed. Indeed, regarding the stability of this inverse problem, 
Alessandrini proved that, assuming n > 3 and a-priori bounds on 7 of the 
form 

Tl 

(1.4) ll7llffo(n) < forsomes > --f 2, 

7 depends continuously on A.^ with a modulus of continuity of logarithmic type. 
We also refer to [Alj . |A2) . which improve the result in m for conductivities 
7 e IT2'°“(B). We refer to |R-B-R| . |B-E-R| and [Liu] for the two-dimensional 
case, where logarithmic type stability estimates have been established too. The 
common logarithmic type of stability cannot be avoided |A31 IMaj . However, the 
ill-posed nature of this problem can be modified to be conditionally well-posed 
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by restricting the conductivity to certain function subspaces. Well-posedness is 
here expressed by Lipschitz stability. 

A first result of this kind was established by Alessandrini and Vessella |A-V) . 
to which we refer, together with m, for an in-depth description and analysis. 
The result of |A-V) was extended to different types of problems, for example, 
in |Be-dH-Q] , |Be-dH-Q-S| for the Schrbdinger and the Helmholtz equations, 
respectively, in [Be-Frj for the inverse conductivity problem with complex con¬ 
ductivity, and in [Be-Fr-V] , [Be-Fr-Mo-Ro-Ve] for the determination of the Lame 
parameters in the elastostatic problem. All of these papers have in common the 
stable determination of coefficients that are piecewise constant on a given, that 
is, fixed domain partition. This partition needs to satisfy certain geometric con¬ 
ditions. One can view the domain partition as prior information which needs to 
be obtained by other inverse methods. 

Let us emphasize that the main effort in these papers resides in achieving a 
constructive evaluation of the Lipschitz constant. This goal requires the con¬ 
struction and evaluation of ad hoc singular solutions and the use of quantitative 
estimates of unique continuation, and both steps may be somewhat sophisticated 
due to the presence of jumps in the coefficients. These are common themes of 
the above mentioned papers, each of which however presents its own specificity 
and difficulty. Also some variations from the route first outlined in |A-V] have 
been followed, let us mention for instance E e-Fr-V| . |Be-dH-Q-S| , by taking 
advantage of a general functional analytic framework developed in [Bac-Vj . 

We wish to recall, here, that the uniqueness result obtained by Druskin |D2| was 
in the context of piecewise constant conductivities too. In the present paper, 
we consider conductivities that are piecewise linear instead. We note that we 
can adapt our analysis to the case of piecewise linear resistivities. 

In dimension n > 3 - which we consider in geophysics - uniqueness has been 
established by Haberman and Tataru for conductivities in iHa^T] and more 
recently for Lipschitz conductivities in |Ca-R] . both assuming full boundary 
data. The original uniqueness result by Sylvester and Uhlmann |Sy-U| required 
the conductivity to be C°°. For the two-dimensional case we refer to |Bro-U] 
and the breakthrough paper |As-P] where uniqueness has been proven for con¬ 
ductivities that are merely L°°. 

The class of conductivities considered in this paper consists of piecewise linear 
functions on a given domain partition, which are possibly discontinuous at the 
interfaces of this partition. The Lipschitz stability estimate we provide requires a 
direct proof. (Indeed, the uniqueness result of |Ca-R) (n > 3) does not apply; in 
fact, in the case of partial data, the result of |As-P] does not apply either). This 
estimate is given in terms of the local DtoN map. With a slight modification, 
our arguments can apply when the local NtoD map is available instead, see for 
instance the discussion in [A-Glj . 

With a Lipschitz stability estimate at hand, we can apply certain iterative meth¬ 
ods for reconstruction within a subspace of piecewise linear functions with a 
starting model at a distance less than the radius of convergence to the unique 
solution |dH-Q-S| , [L-Sj . This radius is roughly inversely proportional to the sta¬ 
bility constant appearing in the estimate. More importantly, we can iteratively 
construct the best piecewise linear approximation for a given domain partition. 
Since the stability constant will grow at least exponentially with the number of 
subdomains in the partition m. the radius of convergence shrinks accordingly. 
One can expect accurate piecewise linear approximations with relatively few 
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subdomains to describe the subsurface, noting that the domain partition need 
not be uniform and may show a local refinement, and hence our result provides 
the necessary insight for developing a practical approach with relatively minor 
prior information. Whether we can recover, also, an unknown domain partition 
(such as one of tetrahedral type) is a current subject of research. 

As we mentioned earlier, the application we have in mind here is the DC acqui¬ 
sition and method, which were introduced by Schlumberger in 1920 [S]. Initial 
DC deep resistivity studies of Earth’s crust were carried out as early as in 1932 
[SS] . Many studies have followed. We mention, in particular, the experiments 
and results by Constable, McElhinny and McFadden |Co-McE-McF) carried out 
in central Australia in 1984. For a general description and the history of the DC 
method (and the closely related induced electrical polarization (IP) method) we 
refer to the textbooks by Koefoed [K2] . Zhdanov and Keller [ZKj . and Kauf¬ 
man and Anderson [Ko-An] : for a concise tutorial and review, see Ward [Wj . 
For a finite-element method and solver for and computational studies of the 
DC method, see Li and Spitzer lESl- Here, we consider isotropic conductiv¬ 
ities (and therefore resistivities); however. Earth’s materials can certainly be 
anisotropic, which was already recognized by Mallet and Doll |M-D| . We refer 
to [AS], [ASI], |As-L-P| . HI, |F-K-R| . [GS], [SE], 10 and [IAMT] ! for results 
concerning the anisotropic case. 

Through recent decades, electromagnetic methods have been widely used in 
geothermal prospecting |Br-Ma-V-F-Mo-E| . Amongst different geophysical ex¬ 
ploration methods, in geothermal prospecting, resistivity methods have been 
demonstrated to be the most effective. The reason is that the electrical re¬ 
sistivity of rocks is controlled by important geothermal parameters including 
temperature, fluid type and salinity, porosity, permeable pathways, fracture 
zones and faults (structural), the composition of the rocks, and the presence of 
alteration minerals. In this context, we mention the work of Hersir, Bjornsson 
and Eysteinsson |He-Bj-E] and, more recently, of Flovenz et al. [FI] . 

We briefly mention how the acquisition - essentially probing the NtoD map 
- is carried out (see, for example, [Balj . [Ba2j l. The original acquisition was 
designed for two-dimensional configurations (n = 2). The Schlumberger ar¬ 
ray consists of four collinear electrodes. The outer two electrodes are current 
(boundary source) electrodes and the inner two electrodes are the potential (re¬ 
ceiver) electrodes. The potential electrodes are installed at the center of the 
electrode array with a small separation. The current electrodes are gradually 
increased to a greater separation during the survey - while the potential elec¬ 
trodes remain in the same position until the observed voltage becomes too small 
to measure - for the current to probe deeper into the earth. Indeed, the depth 
resolution of the DC method is sensitive to the separation between current 
electrodes [QH. There is also the (crossed) square-array acquisition which is 
designed to be more sensitive to anisotropy than the Schlumberger array [tFW| . 
0 . 

There are different types of electrode configuration that are commonly used. 
In two-dimensional configurations, the dipole-dipole array is widely being used 
because of the low electromagnetic coupling between the current and potential 
circuits. In three-dimensional configurations, the pole-pole electrode configura¬ 
tion is commonly used. (In practice, the ideal pole-pole array, with only one 
current and one potential electrode does not exist. To approximate the pole- 
pole array, the second current and potential electrodes must be placed at a 
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large distance). For convenience the electrodes are arranged in a square grid 
with the same unit electrode spacing in orthogonal (coordinate) directions. (We 
mention the E-SCAN method |L-0) . [E-0) 1. It can be very time-consuming to 
make such a large number of measurements. To reduce the number of measure¬ 
ments required without seriously degrading the resolution, “cross-diagonal sur¬ 
vey” method was introduced; here, the potential measurements are only made 
at the electrodes along two orthogonal directions and the 45 degrees diagonal 
lines passing through the current electrode (extracted from Loke’s tutorial: 2-D 
and 3-D electrical imaging surveys, IE]). 

The inverse problem pertaining to resistivity interpretation was reported as 
early as the 1930s (e.g. Slichter, 1933; Stevenson, 1934; Ejen, 1938; Pekeris, 
1940 [Ei]). Slichter jSl] published a method of interpretation of resistivity data 
over a planarly layered earth using Hankel’s Fourier-Bessel inversion formula. 
It gives a unique solution if the resistivity is a continuous function of electrode 
spacings. A substantial number of papers have been written on approaches 
based on partial boundary data “fitting” or optimization to estimate the re¬ 
sistivity, without knowledge of uniqueness or convergence. Narayan, Dusseault 
and Nobes [Na-D-No] give an extensive overview. In the context of data fitting, 
Parker [P] indicates and illustrates in planarly layered models the Ill-posedness 
of the IBVP. Interestingly, various studies and implementations have resorted to 
“blocky” (and pseudo-layered) representations of resistivity i |Lo-A-n| . [EQ], 
[Au-VC] ') and, hence, fit the class of functions for which Lipschitz stability esti¬ 
mates have been obtained. Finally, we mention the complementary frequency- 
dependent transient electromagnetic (TEM), magnetotelluric (MT) and electro- 
seismic methods. The hybrid inverse problem of electroseismic conversion was 
analyzed by Chen and De Hoop [Ch-dH] . The further analysis of TEM/MT 
m is a subject of current research. 

In recent years, there has been a renewed and growing interest in the application 
of electrostatic and diffuse electromagnetic inverse boundary value problems in 
geophysics driven by the idea of combining different probing fields, including 
acousto-elastic waves, to identify the (poro-elastic) rock properties in Earth’s 
interior within a particular geological structure in an integrated fashion. These 
properties certainly will not vary smoothly. We capture the geological structure 
in a domain partition, let the properties be discontinuous across subdomain 
boundaries of geological significance, and approximate the parameters, here 
conductivity in the electrostatic problem, in each subdomain by linear interpo¬ 
lation. (From a rock physics point of view, this interpolation should be obtained 
from a nonlinear upscaling, which is still an active area of research). This ap¬ 
proach, and the generality of these approximations, analyzed in the context of 
conditional well-posedness are the novelty of this paper. 

The outline of the paper is as follows. Our main assumptions and our main 
result iTheorem 12.311 are given in section 2. Section 3 contains the proof of the 
main result, as well as two intermediate results iTheorem 13.21 and Proposition 
13.3[) needed to build the necessary machinery. Theorem 13.21 provides original 
asymptotic estimates for the Green’s function of the conductivity equation, its 
gradient and a mixed derivatives, for conductivities that are linear on each do¬ 
main Dj of a given partition {Dj} of Q. These asymptotic estimates are given at 
the interfaces between the domains Dj, where the conductivity is discontinuous. 
Proposition [331 provides estimates of unique continuation of the solution to the 
conductivity equation for piecewise linear conductivities. Section 4 is devoted 
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to the proof of Theorem 13.21 and Proposition 13.31 The latter is based on the 
argument introduced in [A-V] [proof of Proposition 4.4], therefore only the main 
differences in the two proofs are highlighted in the present paper. 


2 Main Result 

2.1 Notation and definitions 

In several places within this manuscript it will be useful to single out one coor¬ 
dinate direction. To this purpose, the following notations for points x € K" will 
be adopted. For n > 3, a point x € R" will be denoted by a; = (x',Xn), where 
x' G and Xn G R. Moreover, given a point x G R", we will denote with 

Br(x), B'j.(x) the open balls in R", R"“^ respectively centred at x with radius r 
and by Qr(x) the cylinder 

Qr{x) = B'^{x') X {Xn -r,Xn+ t). 

We will also denote 


R!^ = {ix',Xn) G R’^|x„ > 0}; R” = {{x',Xn) G R"|a;„ < 0}; 

B+ = Bril R!]:; B- = BrllWl; 

g+ = g,nR"; Q; = QrnRi. 

In the sequel, we will make a repeated use of quantitative notions of smoothness 
for the boundaries of various domains. Let us introduce the following notation 
and definitions. 

Definition 2.1. Let LI be a domain in R”. We say that a portion S of dLl 
is of Lipschitz class with constants tq^L if for any P G dT, there exists a rigid 
transformation of R" under which we have P = 0 and 

LI n Qro — G Qvq • ^ 

where ip is a Lipschitz function on B'r^ satisfying 

(p(0)=0; \\ip\\co,i(^B'^^) < LrQ. 

It is understood that dLl is of Lipschitz class with constants ro, L as a special 
case of S, with E = dLl. 

Definition 2.2. Let LI be a domain in R". We say that a portion E of dLl is 
a flat portion of size rg if for any P G T, there exists a rigid transformation of 
R" under which we have P = 0 and 

E n Qro/3 “ L. Qro/sl^n — 0} 

LI n Qtq/S — L Qro/sl^n ^ 0} 

(2.1) (R" \ n) n Qro/S = {x L Qro/flXn < Ll}, 

Let us rigorously define the local D-N map. 
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Definition 2.3. Let Q be a domain in K" with Lipschitz boundary dLl and S 
an open non-empty (flat) open portion of dil. Let us introduce the subspace of 

H^dn) 

( 2 . 2 ) i?i(E) = {/ e H^dn) I supp/ c s}. 

_ _1 

and its dual Hco^ (S). Assume that 7 S L°°{Ll) satisfies 

(2.3) < 7 ( 0 :) < A, for almost every x € ft, 

then the local Dirichlet-to-Neumann map associated to 7 and S is the operator 

(2.4) : i7i(E) ^ H~^ (S) 
defined by 

(2.5) < g, r] > = / '^{x)Vu{x) ■ V4>{x) dx, 

J n 

for any g, t} ^ where u G is the weak solution to 

( div{’y(x)\/u(x)) = 0, in LI, 

\ u = g, on dfl, 

and (j) € H^{Ll) is any function such that 4>\dn = g in the trace sense. Here we 

1 _1 

denote by < ■, ■ > the Lp‘{dLf)-pairing between and its dual Hco^lYf). 

Note that, by ( 1 ^ . it is easily verified that is selfadjoint. We will denote 
by II • II* the norm on the Banach space of bounded linear operators between 
i?i(E) andiJ’^S). 

- 

Remark 2.1. Note that the space RQ^g(I]) (\LiMf . Chapter 1) is the closure 

i I 

of HfofCf) in H"i{dLl), therefore the local DN map could be equivalently given 

- —- - 

by replacing in definition \2.S\ the spaces Hfo{'Lf), Hco^(Yf) with Hqq(T,) and its 
_ _1 

dual Rqq^ (S) respectively and by continuing to use the notation < •, • > for the 

1 _1 

Lf{dLl)-pairing between Hqq{'S,) and (S). 

2.2 Assumptions 

2.2.1 Assumptions about the domain Ll 

1. We assume that 11 is a domain in R" satisfying 


(2.6) |fl| < Nr^, 

where |fl| denotes the Lebesgue measure of Ll. 

2. We fix an open non-empty subset S of dLl (where the measurements in 
terms of the local D-N map are taken). 
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N 

n = u 4, 

i=i 

where Dj, j = 1,... ,N are known open sets of ffi.", satisfying the condi¬ 
tions below. 

(a) Dj, j = 1,..., N are connected and pairwise nonoverlapping polyhe¬ 
drons. 

(b) dDj, j = 1,... ,N are of Lipschitz class with constants tq, L. 

(c) There exists one region, say Di, such that dDinY, contains a flat por¬ 
tion Si of size ro and for every i G {2,..., N} there exists ji,... ,jK G 
{!,... ,7V} such that 

(2.7) Dj^ = Di, = Di. 

In addition we assume that, for every k = 1,..., iV, dDj^ fl dDj^_^ 
contains a flat portion Efc of size rp (here we agree that Dj^ = R"\n), 
such that 

S/c C n, for every k = 2,..., K, 

and, for every k = 1,..., K, there exists Pk G Sj, and a rigid trans¬ 
formation of coordinates under which we have Pk = 0 and 

Sfc l~l Qro/S ~ ^ Qro/sl^^n “ 0 } 

^jk ^ Qro/3 ^ Qro/sl^^n ^ 0} 

( 2 . 8 ) ^jk-l fl Qro/3 “ S Qro/sl^n ^ 0 }) 

2.2.2 A-priori information on the conductivity 7 

We will consider a conductivity function 7 of type 

N 

(2.9a) 7 (a:) = ^ lflx)xD, {x), xGO., 

i=i 

(2.9b) Ijix) = aj + Aj ■ x, 

where aj G R, Aj G R" and Dj, j = 1,... ,N are the given subdomains intro¬ 
duced in section [ 2 . 2. II We also assume that 

(2.10) < 7 j < A, a.einfl, for any j = 1,... n, 

for some positive constant A. 

Definition 2.4. Let N, tq, L, M, A be given positive numbers with A G N. 
We will refer to this set of numbers, along with the space dimension n, as to the 
a-priori data. 
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Remark 2 . 2 . Observe that the class of functions of the form (|2.9al) - (I2.9bl) is 
a finite dimensional linear space. The L°° - norm || 7 ||L=o(n) is equivalent to the 
norm 

IIItIII = maxj=i,...,Ar{|aj| + |Aj|} 

modulo constants which only depend on the a-priori data. 

From now on for simplicity we will write 

A, = i = l,2. 

Theorem 2.3. Let LI, Dj, j = 1.,..., N and H be a domain, N subdomains of 
n and a portion of dLl as in section V2.2.1\ resvectivelu. Let i = 1,2 be two 
conductivities satisfying ( 12 . 101 ) and of type 


N 

( 2 . 11 ) = J2ij"\x)xD,{x), xen, 

1=1 

where 





■X, 


with € K and € R", then we have 


(2.12) <C||Ai-A2||., 

where C is a positive constant that depends on the a-priori data only. 

Remark 2.4. In this paper we are assuming that the conductivity "f is a piece- 
wise linear function. However, the case where the resistivity function p = 7 “^ 
is piecewise linear, could be treated equally well. 

3 Proof of the main result 

The proof of our main result (theorem 12.3|) is based on an argument that com¬ 
bines asymptotic type of estimates for the Green’s function of the operator 

(3.13) L = div ( 7 (a:)V) in H, 

ftheorem 13.2|) . with 7 satisfying (I2.9al) - (I2.10|) . together with a result of unique 
continuation (proposition [331) for solutions to 

Lu = 0, in n. 

Our idea in estimating 7 ^^^ — q(2) exploits, on one hand, an estimate from below 
of the the blow up of some singular solutions (which we will introduce below) 
Su and some of its derivatives if 7 ^^^ — 7 ^^^ is large at some point. On the other 
hand, we will use estimates of propagation of smallness to show that Su needs 
to be small if Ai — A 2 is small. We will give the precise formulation of these 
results in what follows. 
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3.1 Singular solutions 

We will start with some general considerations about the Green’s function 
G(x, y) associated to the operator (Id.ldL where 7 is merely a measurable matrix 
valued function satisfying the ellipticity condition (ESI). 

3.1.1 Green’s function 

If L is the operator given in (I3.13p . then for every y G ft, the Green’s function 
G{-,y) is the weak solution to the Dirichlet problem 


/ 01 ', / div( 7 VG(-,y)) =-d(-- y) , in 12 , 

1 G(-,y) = 0 , on 512 , 

where 6{- — y) is the Dirac measure at y. We recall that G satisfies the properties 
l |T.it-St-Wp 


(3.2) G{x,y) = G{y,x) foreveryx, y e 12, x ^ y, 


(3.3) 0 < G(x, y) < G|x — yp for every x, y e 12, x^y, 

where G > 0 is a constant depending on A and n only. Moreover, the following 
result holds true. 

Proposition 3.1. For any y S 12 and every r > 0 we have that 

(3.4) [ |VG(-,y)p <Gr2-" 

Jn\Br(y) 

where G > 0 depends on A and n only. 

Proof. The proof can be obtained by combining Caccioppoli inequality with 
(13.311 1 [A-V) . Proposition 3.1). □ 

3.1.2 The Su singular solutions 

Let 7 ^*^ z = 1, 2 be two measurable functions satisfying the ellipticity condition 
(E^l) and let Gi{x,y) be the Green’s functions associated to the operators 

(3.5) Li = div ^ 7 *'*^(x)V^ in 12, z = l,2. 

Let U be an open subset of 12 and W = 12 \ G. For any y, z G W we define 

(3.6) Su{y,z)= ( {'y'^'^\x) - j^‘^'>{x))Va;Gi{x,y) ■Va:G 2 {z,x)dx. 

Ju 
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We recall that (see [A-Vj 'l for every y,z £ W we have that Su{-, z), Suiy, ■) £ 
Hlo^{W) are weak solutions to 


(3.7) 

div (j^^'>{-)VSui-,z)^ = 0, 

in W 

(3.8) 

div (j^‘^H-)ySuiy,-)^ = 0 , 

in W. 

It is expected that Suiy, z:) blows up as y, z approach simultaneously one point 
oidU. 

We will denote with 

(3.9) 

r(a;,y)-. . |a; y\ 

[n - 2)0Jn 



the fundamental solution of the Laplace operator (here w„/n denotes the volume 
of the unit ball in ffi.”). li Di, i = 1,... ,N are the domains introduced in section 
12.2.11 and L is the operator given by (13.1311 . we will give asymptotic estimates 
for the Green’s function of L, with respect to (13.9|) at the interfaces between the 
domains Di, i = 1,... N. These estimates are given below. In what follows let 
G be the Green’s function associated to the operator L in 17. 

3.1.3 Asymptotics at interfaces 

Theorem 3.2. Let Qi+i be a point such that Qi+i £ Sio (7^i+i) H ^i+i with 
ZG{l,...,iV — 1} . There exist constants l3,0,0</3 <1,0 <0<1 and C > 0 
depending on the a priori data only such that following inequalities hold true for 
every x £ (Qi+i) H and every y = Qi+i — re„, where r £ (0, ■^) 


(3.10) 

(3.11) 

(3.12) 


G{x,y) - 


V,G(i,y)- 


Iji iQi+i) +7ii+i(Q/+i) 
2 


r(a;,y) 


< G\x - y\ 


3 — n 


VyVxG(a;,y) - 


IniQi+i) + 7ii+i(<3i+i) 
2 


Va;T{x,y) 


<G\x-yf+^-^, 


-VyVj:T{x,y) 


< G\x - y\ 


9—n 


3.2 Quantitative unique continuation 

We recall that up to a rigid transformation of coordinates we can assume that 

Pi = 0 ; (R" \ 17) n Bro = {{x', Xn) £ Pro I Xn < T[x')}, 

where (/? is a Lipschitz function such that 

(^(0) = 0 and ||¥’||co.i(b;j < Lr^. 

Denoting by 


Dq = lx £ 


\ 17) n Pr 


< g^’o, * = l,...,n- 1, 


xo 
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it turns out that the augmented domain ^Iq = flU Dq is of Lipschitz class with 
constants ^ and L, where L depends on L only. We consider the operator Li 
given by (13.511 and extend 7 ^*^ to 7 ^*^ on LIq, by setting 7 ^*^|_Do = 1, for i = 1,2. 
We denote by Gt the Green function associated to Li = div( 7 (*^(a;)V-) in LIq, 
for i = 1,2. For any number r € (O, |ro) we also denote 

{Do)r = {x G Dq I dist{x, n) > r} . 

Let K G 1,..., N he the subdomain of ft such that 

(3.13) E = 117 ( 1 ) - 7^^^llL“(n) = - 1^^^\\l^{Dk)- 

and recall that there exist ji, - ■ ■ ,jK G 1,... ,N such that 

Dj^ E\ ,... Ej^ , 

with ,... satisfying assumption 4(ci). For simplicity, let us rearrange 
the indices of these subdomains so that the above mentioned chain is simply 
denoted by Di,..., Ek, K < N. We also denote 


(3.14) yVk = [}E„ Uk=Eo\Wk, forfc = l,...,iF 

(3.15) ~Su,{y,z)= j ( 7 (')- 7 ^"))VGi(-,y)-VG 2 (-,z), for A: = 1,..., X. 

JUk 

We introduce for any number 6 > 0 as in [A-V] . the concave non decreasing 
function ujb(t), defined on (0,+oo), 


UJb{t) = 


2''e-2|log^|-^ tG(0,e-2) 


0-2 


t G [e ^, + 00 ) 


and denote 

(3.16) 

The following parameters will also be introduced 


Wf, — UJ, U)), — Wb o Wf, 


P = arctan —, 

1j 

Pi = Ai sin ,01, 


Pi = arctan 

1 — sin Pi 
1 + sin Pi 



Ai 


1 + sin pi 


\m = aXm-i, Pm = apm- 1 , foreverym>2, 

dra — XjYi Pm: Xfl ^ 1 . 

For k = 1,... ,K and a fixed point y G Sfe+i, denote 


(3.17) Wm{y) = y - Xmv{y): for every TO > 1 , 

where ^{y) is the exterior unit normal to dEk- The following estimate for 
(y, z) holds true, for fc = 1,..., AT. 
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Proposition 3.3. (Plstimates of unique continuation^ If, for a positive 
number So, we have 


(3.18) Suk{y,z) < rp "-£ 0 , for every {y,z) e {Do)ra. X iDo)rf., 


£0 \\ 


then the following inequality holds true for every r G ( 0 ,di] 

(3.19) 

Su, {w-^iQk+i),wj,{Qk+i)) < r-^+^C^{E + eo) j j 

(3.20) 

dyA.Su, iwj,iQk+i),y^hiQk+i))\ < rAC^E + eo) 


{i/cf 


(l/C)'- 


for any i,j = where Qk+i G Sfc+i fl ilm(Pfe+i), h{r) = min{m G 

N I dm < r}, w;,(^)((3fe+i) = Qk+i - A^(r)i/(Qfe+i), Am has been introduced 
above, v is the exterior unit normal to dDk and Ci,C 2 > 0 depend on the 
a-priori data only. 


3.3 Lipschitz stability 

Proof of Theorem \2.tA Let Dk be the subdomain of V, satisfying (I3.13|) and 
let Di,...Dk be the chain of domains satisfying assumption 4(d). For any 
k = 1,... ,K we will denote by Dxf and 9^/ the n — 1 dimensional vector of 
the tangential partial derivatives of a function / on and the normal partial 
derivative of / on respectively. Let us also simplify our notation by replacing 
with Ai, for i = 1, 2. We will also denote 

£0 = l|Ai - A 2 II*, Si = - 7 ^^^ 1 1 LOO (Wi)- 

We start our argument by estimating di. By |A-V) we obtain the following 
estimate on the first interface Ei (this can also be obtained as a straightforward 
consequence of |A-G1] [Theorem 2.3] for 7 *^^) — continuous near Ei) 


(3.21) 117^1 )-7 


( 2)1 




< C(eo + E) ( 


£0 


£0 + if 


where ro > 0 is the constant introduced in subsection l2.1l and C > 0 is a constant 
depending on the a-priori data only. Let us denote by 


(3.22) ai-f/3i •a;= ( 7 ^^- 7 f^)(a:), 

and by {ej}j=i^,,,^n-i a family of n — 1 orthonormal vectors starting at Pi, 
defining the hyperplane containing the flat part of Ei. By computing 7 ^^^ — 7 p^ 
on the points Pi, j = 1 ,..., n— 1 , taking their differences and applying 

(13.211) . we obtain 
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(3.23) 


|ai+/3i-Pi| < C{eo + E) 


C 


(3.24) 


to 


,(2) j 

f £0 

1/C 1 

K,£o + E 

,(2) j 

( £0 

1/C 1 

\£q E 


for j = 1 , where C > 0 is a constant depending on the a-priori data 

only. To estimate fii along the remaining direction v and therefore ai, the 
normal derivative 


needs to be estimated. We recall that for every y,z € Dq we have 


((Ai-A 2 )Gi(., 2 /),G 2 (.,z)) 

(3.25) 

and 


/ (7^')-7^"))(-)VGi(-,y)-VG2(-,z) 
JQ, 

Suo{y,z), 


((Ai-A2)9,„Gi(-,y),a,„G2(-,z)) 

(3.26) 

From (13.2511 we obtain 


/ ( 7 ^'^ - 7 ^"^)(-) 9 y„ VGi(-, y) ■ VG 2 (-, z) 

JQ, 

dy„dzr,Suo{y,z). 


(3.27) 


\Suo{y,z)\ 


< Geo?'o~"; for every G (£>o)r'o/3, 


where G depends on A, L, X and n. Let po = where G is the constant 
introduced in Theorem 13.21 let r G ( 6 ,^ 2 ) and denote 


w = Pi + av, 


where CT = a^ ^Ai, 


then 


(3.28) dyndzr.Suoiwy'w) = h{w) + hiw), 

where 

Ii{w)= [ (7''^-7^'^)(-)5y„VGi(-,u;)-a,„VG2(-,ri;), 

J BpQ (Pi)nDi 

l2{w) = [ (7(i)-7(2))(-)ay„VGi(-,u;)-a,„VG2(-,a;) 

Jn\{Bp^{Pi)nDi) 

and (see [A-V ] ) 
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(3.29) |/ 2 (zi;)| < CEp^\ 

where C depends on A and n only. We have 




> 


' BpQ (Pi)nDi 


- Pi)ndy„VGi{-, w) ■ d,^VG 2 {-,w) 


BpQ (Pi)nDi 


' BpQ (Pi)nDi 


- 7^"))(Pl)) • (x - Pi)'P,„VGi(-,ra)| VG2(-, «;)| 

K^W-^f))(Pi)P^^VGi(.,«;)| |a.„VG2(-,u;)|. 


and, by Theorem 13.21 this leads to 


|/i(w)| -ji^^)\Gi [ \dy„V^T{x,w)f\xr, 

J BpQ (Pi)nPi 


X —w\ , 


-C 2 E / |5,„V,r(a;,^z;)| 

J Bpg{Pl)nDl Pq 

f |a; — 


-G.E 


, s 2/3 

' BpgiPk)nDi Pq 

„( 1 ) _ .,( 2 ) 


Ppo (Pi)nDi 


|i^r(7r -7101 \i^-Piy\ |9,„VGi(.,«;)| |5.„VG2(-,rr;)| 


(3.30) - / 1(7^^ -7f^)(i^i)l |5,„VGi(-,«;)| \d.^VG2{;w)\, 

J Ppp (Pi)nPi 

where Ci, 6 ^ 2 , Ca are constants that depends on M, A and n only. Therefore, by 
combining p.30l) together with (13.281) and (13.291) . we obtain 


|/i(w)| > \du{li^ - li'’)\Ci I \x-w\ 

BpQ (Pi)nDi 

C 2 E f 


which leads to 


\x — w 


l-2n+/3 


Po Bpg{Pl)nDl 

CqE f 

~ ^ . I"" 

Po J Bpg(PP)^^D^ 

— £oG4 / \x — W 

Bpg (Pl)nPl 

- eoGs [ 

J PpQ (Pi)nPi 


,|2 —2n+/3 


l-2n 


(3.31) 19^7^^ - 7f < |A(3«)| + Gieoa-^ + G 3 


1—n+/3 


Po 
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where 


(3.32) |/iH| < \dy^d.Mw,w)\ + CEp^-. 

Thus by combining the last two inequalities we get 


-n+/3 


191,(71^^-7f^)|cr^ ” < \dy^d^^Suo{w,w)\+CE'^-^ 

Po 


+ + C 3 - 


(3.33) 

and by recalling that by Proposition 13.31 we have 
dy, d,, Su, {w, w) < {E + eo) 


we obtain 
(3.34) 

|5.(7r^-7f^)l < ( C^^''\E + e^) 


jl-n+P 

Po 


£0 \ 


( 1 /C)'*M 


E + eo 


£0 \ 


(l/C)hM 


E + £0 


cr^\ ~ 

+ Cieo + CE^\+C3 — 

Po J P 


We need to estimate and in terms of r, where C > 1. It turns out 

that 


2 / di \ Toiy 


Cn < —) 

l.h 


< 


(3.35) (-) 
therefore for any r G (0, (^ 2 ) 

(3.36) 

|5.(7r^-7f^)l < (^) C(i?+eo) 

which leads to 


1 \2 / r \-is^ 


(^)(i 


£0 


E + Eq 




+ 


di J I \di 


£ 0 , 


(3.37) 
and thus 


l'9i-(7f ^ - 7f^)l < C{eo + E) ( 


( 2 ) / £0 


(3.38) 


1^1 .^<^(.0 + ^;) 
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By combining (13.381) together with (I3.22L (13.231) and (I3.24L we get 


(3.39) Kl, j/3,1 < C(eo + E) 
and by (13.391) 

(3.40) < C(eo + E) " , 

where C > 0 is a constant that depends on the a-priori data only. By proceeding 
by induction on I in order to estimate 7 ^^^^ — for I = 1,... ,K, we replace 
(13.251) and (13.261) by 


((Ai-A2)Gi(.,y),G2(-,z)) 

(3.41) =Su,_Ay,^)+ [ (7^')-7^"P(-)VGi(-,y)-VG2(-,z) 

JWi-i 


and 


((Ai-A2)a,„Gi(-,y),a.„G2(-,0)) 

{3A2) =dy^d.Ju,_Ay,z)+ [ (P')-7^"P(-)a,„VGi(-,y)-a.„VG2(-,z) 

JWi-i 

respectively. By noticing that (13.411) leads to (see [A-V) l 

(1/5- 2^)1 ^ 11^1 (‘5 y)11^1/2^^pG2(*, 2;)11^1/2^^^ 

(3.43) < C(eQ+5i-i)rl~'^, for every ?/,z G (i:»o)ro/ 3 , 

where G depends on A, L, A, n and by repeating the same argument applied 
for the special case I = 1 and observing that 

5i < Si-i + ||7p^ - 7p^l|L-(Di)> 

we obtain (see |A-V) 1 

< 6l.^ + G(eo + Sl.^ + E) (_££±^^) ^ , 

in particular for I = K 

Sk <Sk-i+ C{so + Sk-i + E) ^ , 

which leads to 
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therefore 


Ih'" - 7‘">lll.(0, < C(e + E) (uf' (^)) ° ■ 


(3,44) B<C(.. + £)(.f 

Assuming that E > Sge^ (if this is not the case then the theorem is proven) we 
obtain 


which leads to 



therefore 


E < 


1 


LO 1 


(^) 


£ 0 ) 


where here, with a slight abuse of notation, cu\ ' denotes the inverse function 

ZT 

of a;\ . 

c 

□ 


4 


Proofs of Theorems 


3.2 


3.3 


4.1 Asymptotic estimates 

Theorem 4.1. Let r > 0 be a fixed number. Let U £ H^{Qr) be a solution to 


(4.45) 


div{b{x)VU) = 0 , 


where 


b{x) 


b++B+-x, X £ (5+, 
b~ + B~ ■ X, X £ Q~ , 


where b'^,b~ £ M.,B'^,B~ £ R" and 0 < b~^ < b{x) < b. Then, there exist 
positive constants 0 <a'<l,C'>0 depending on b,r and n only, such that for 
any p < § and for any x £ Qr- 2 p, the following estimate holds 


l|Vt/|U~(Q,(x)) +p“'|Vt7U, +P“'|Vt7U, 

,Qp{x)r\Qr 


- pi.^/2 ll^ll^"(Q2pW) 
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Proof. For the proof we refer to |Lad-Ur[ Theorem 16.2, Chap.3], where the 
authors obtained piecewise C^’°‘ estimates for solutions to linear second order 
elliptic equations with piecewise Holder continuous coefficients and discon¬ 
tinuity interfaces (see also flTVn] |in^ for more recent results under weaker 
regularity hypothesis) . □ 

Proof of Theorern, \S.2\. We fix^ G {l,...,iV — 1}. With no loss of generality we 
may assume that Qi+i = 0, a; = 1 and ai+i > 0. We will denote a; = a~ and 
o-i+i = 

For any x = {x', Xn) we denote x* = {x', —Xn) and we have that a fundamental 
solution of the operator diva;(l + (a’*’ — l)x^V 2 ,) has the following explicit form 

[ y) + a4-°(a++i) r(a^, y*) , if ?/„ > o 

(4.47) H{x,y) = l ■^^T{x,y) , if < 0 

[ r(a;,j/)-h i^r(a;,y*) , ifa;„,j/„< 0 . 

We then define 

(4.48) R{x, y) = G{x, y) - H(x, y) = G{x, y) - ^ r(x, y) . 

We observe that R in (14.481) satisfies 

(4.m / =-divx(( 7 (-)- 7 o(-))Va;id(-,y)) , in H , 

^ \ R{;y) = -Hi;y) , on dQ . 

By the representation formula over H we have that R in (14.4811 satisfies 

(4.50) R{x, y) = - f (7(C) - 7o(C))VcH(C, v) ' VcG(C, x)dC 

JQ 

(4.51) + [ 7(C)a.G(C,x)i7(C,j/)da(C). 

Jdn 

We hrst treat the boundary term on the right hand side of the above equation. 
We have that 


(4.52) 

(4.53) 

(4.54) 


G(-, x) H{-, enyn) da 


< 7 

< 7 l|G(->a:)||^i(n\B,.^^2(a::)) ■ 


Hence, we deduce that 


(4.55) 



G(-,x) i7(-,e„z/„) da 


< G 


where G > 0 is a constant depending on the a priori data only 
Being 7(^) — 70(C) of Lipschitz class, we observe that the estimate 
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(4.56) I /(7(C)-7o(C))Vci?(C,e„?/„)-VcG(C,x)dCI 

Jn 

can be achieved along the lines of the proof of Claim 4.3 in [A-V] . 

Combining (14.551) and (14.561) we get 

(4.57) \Rix,enyn)\ < Gi|a; - , 

when X G and yn G {-ro,0). _ 

We now focus on the estimate for Va:i?(a;, e„y„). Again arguing as in [A-VI 
Claim 4.3], we fix a; G and yn G (— ro/4,0) and let us denote 

(4.58) Q = B'f^/^{x') X (^Xn,Xn + ■ 

where h = \x — y\ . We observe that Q C Moreover, we have that 

2 

Q C Qh{x) and X G dQ. 

By p.3l) . Theorem 14.11 and explicit computation on the behaviour of H{x, y) we 
get 


(4.59) |V,G(-,e„2/„)U^Q , |V,iJ(-, e„2/„)U-,Q < Ch-^'+^-^ 

where G > 0 is a constant depending on the a priori data only. 
Hence by (14.481) and (14.591) we get 


(4.60) |V,i?(-, enyn)U',Q < G/i-“'+i-" . 

where G > 0 is a constant depending on the a priori data only. 
We recall the following interpolation inequality 


(4.61) ||Va;i?(-,e„y„)||ioo(Q) <G||i?(-,e 


7/ iir'G+“' 


v,i?(- 


77 ')|G1+" 
T-yn)\a' ,Q 


where G > 0 is a constant depending on the a priori data only. 
By the above estimate and (14.571) we obtain 

(4.62) \V^R{x,y)\<Ch^+^-^ 


where /? = ^ . 

Finally, we study the behaviour of S/yVxR{x,y). Let us define the cylinder 
Q = B'hiO) X {yn — j,yn)- As before we have that Q C Q~ C Qh{y). In 
particular, we have that x ^ Qh{y)- 

4 

Let k be an integer such that k G {1,... ,n}. We observe that dxt,T{x, •) and 
dx^G{x, ■) are solutions to 

Ay((9,cj,r(a;, •)) = 0 in Qh{y) , 

Ah7y{-^{-)Vydx^G{x,-)) =Q in Qh(7/) 


(4.63) 

(4.64) 
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respectively. 

Again by applying Theorem 14.11 we have that 

(4.65) \yyd,Mxr)\^, Q ■ 

4 

where C > 0 is a constant depending on the a priori data only. 

We now fix rj € Qh{y) and we notice that f] ^ Qj!_(a:)). By Theorem 14.11 we 

4 16 

have that 


(4.66) 11 V,G(-,ry)|Uo.(Q , ||G(-,r;)|U«=(Q , < Ch^-^ 

^ TE 

where (7 > 0 is a constant depending on the a priori data only. 

Combining (I4.65|l and (14.6611 we have 

(4.67) |V,5,,G(cr,.)U,^<Gh-“'-" , 

where G > 0 is a constant depending on the a priori data only. By explicit 
computations we infer that 

(4.68) lVyd,,r(x,-)i.^Q<Ch-^'-\ 

where G > 0 is a constant depending on the a priori data only. From (14.671) and 
(14.681) . we have that 

(4.69) |V,a,,i?(a;,.)U,Q<G/i-“'-", 

where G > 0 is a constant depending on the a priori data only. 

Moreover, we observe that by analogous arguments of those discussed above, we 
can infer that 

(4.70) iid,,E(x,-)h^(Q^<Ch^+^-- 

r2 

where /3 = . By the following interpolation inequality 


(4.71) ||V,9,,i?(x,.)||^,.(Q) <G||9,,i?(x,.)||5j'^^|V,a,,i?(x,.)|^ 

and by (14.7011 and (I4.69|) we have that 

(4.72) |V,a,,i?(cr,y)| <Gh-"+^ 

where 9 = , . 


□ 


4.2 Propagation of smallness 

Proof of Theorem AS.S[ By repeating the argument in | A- V) . [proof of Proposi¬ 
tion 4.4] concerning a careful analysis of unique continuation argument across 
K discontinuity interfaces and based on an iterated use of the three spheres 
inequality for elliptic equation, we have that for any y, z G (Qfc-i-i)) 
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/ , ^ / P \ \ (I/C)'*''-) 

(4.73) |5^,(y,z)| <r-+^C'‘W(£; + eo) 

where C > 0 is a constant depending on the a priori data only. Hence (13.191) 
trivially follows from (14.731) . 

We now consider Su^ {y, z) as a function of 2n variables where {y, z) S 
Hence by (|4.73ll we have that 

(4.74) 

..., j/„, zi,..., z„)| < + £o) j ) 

for any a; = (yi,... ,y„,zi,..., z„) S (■u;y(^) (Qfc+i)) x (■u;;,(^) (Q^+i)). 

Now observing that Sukivi, ■ ■ ■, Vn, zi,..., Zn) is a solution in Dk x Dk of the 
elliptic equation 


(4.75) divy(j^(y)VySu,(y,z)) + div^(j^(z)V^SuUy,z)) = 0 

we have that by Schauder interior estimates that for any i,j = 1,... ,n it follows 


Hdy^dzjSu^(yi, . . . ,yn, Zi, . . . , Zn)\\L°°(B (»"fc(r)(Qfc+i))xB pji(r) (™fe(r)(Qfc+l))) 

c 

^ 3 - WSudVl^ . . . ,y„,Zi, . . ■ 

Moreover, we have that being > r, hence it follows r < ■^Ph(r)j which 

in turn leads to 


\\dxidoo^Su,ci^l^ ■ ■ ■ ^ ^'Zn)\\L--{Q p^iw-^^,^{Qk+i))) 

2 

c 

(4-76) < 

where C > 0 is a constant depending on the a priori data only. Noticing that 
(4.77) 

we find that 


log(r/ro) ^ ^ log(^Ao) ^ ^ 


log(a) 


log(a) 


(4.78) 


<1-1 I 4 

ro 


h{r) 


Finally by combining (I4.74L (I4.76P and the above inequality we get the desired 
estimate. 

□ 
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